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ABSTRACT 


The Fermi Large Area Telescope (LAT) light curve repository (LCR) is a publicly available, continu- 
ally updated library of gamma-ray light curves of variable Fermi-LAT sources generated over multiple 
timescales. The Fermi-LAT LCR aims to provide publication-quality light curves binned on timescales 
of 3 days, 7 days, and 30 days for 1525 sources deemed variable in the source catalog of the first 
10 years of Fermi-LAT observations. The repository consists of light curves generated through full 
likelihood analyses that model the sources and the surrounding region, providing fluxes and photon 
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indices for each time bin. The LCR is intended as a resource for the time-domain and multi-messenger 
communities by allowing users to quickly search LAT data to identify correlated variability and flaring 
emission episodes from gamma-ray sources. We describe the sample selection and analysis employed 
by the LCR and provide an overview of the associated data access portal. 


Keywords: Fermi-LAT, gamma-rays, time domain, active galaxies 


1. INTRODUCTION 


Study of variability of astronomical objects has led to 
many discoveries in modern astronomy. Identified as one 
of the central themes of the 2020 Decadal Review, the in- 
creasing realization of time-domain and multi-messenger 
astronomy is now opening an entirely new window on the 
Universe. A high duty cycle and long-term monitoring of 
the gamma-ray sky has made the Fermi Large Area Tele- 
scope (LAT; Atwood et al. 2009) a pivotal tool in the 
study of time-domain and multi-messenger astronomy. 
More than a decade of LAT observations has provided 
the identification and regular monitoring of thousands of 
transient, variable, and steady-state sources (Abdollahi 
et al. 2022a). The long-term monitoring of the gamma- 
ray sky by the LAT played a crucial role in the first as- 
sociation of a high-energy neutrino detected by the Ice- 
Cube neutrino observatory and the flaring blazar TXS 
0506+056 (Aartsen et al. 2018a,b). 

The TXS 0506+056 association provided the first tan- 
talizing clue to the origin of the high-energy cosmic neu- 
trino flux detected by the IceCube neutrino observa- 
tory (Abbasi et al. 2022). Though blazars have long 
been suggested as a possible source of extragalactic 
high-energy neutrinos, constraints on neutrino emission 
from bright gamma-ray blazars have largely disfavored 
steady-state active galactic nuclei (AGN) as the pri- 
mary source of the observed neutrino flux (Aartsen et al. 
2017). The TXS 0506+056 association, however, has 
shown that at least some of the astrophysical neutrinos 
detected by IceCube could be attributed to high fluence 
AGN that undergo periods of intense flaring activity on 
top of a much lower quiescent gamma-ray emission that 
may be undetectable to the LAT. 

The long-term monitoring provided by the LAT has 
also been indispensable to multi-wavelength campaigns 
that aim to study long-term correlated variability in 
AGN. Observations of radio and optical flares have been 
used to examine the location of the gamma-ray emitting 
region, as well as the particle acceleration and radiation 
processes in the relativistic jets in these sources. For 
example, the long-term correlated variability between 
optical and gamma-ray flares has allowed for the de- 
tection of a systematic lag in the optical-to-gamma-ray 
flares in 3C 279 (Hayashida et al. 2012, 2015), the evi- 


dence of quasi-periodic variations in the BL Lac object 
PG 1553+113 (Ackermann et al. 2015), and constraints 
on the rate of orphan optical flares with no gamma-ray 
counterpart (Liodakis et al. 2019). Such studies of cor- 
related multi-wavelength variability are a crucial diag- 
nostic of the physics of relativistic jets. 

Despite the importance of long-term variability stud- 
ies of LAT sources, very few existing resources enable 
the community to easily access the whole-mission data 
to perform correlative analyses. Resources like the Fermi 
Al-Sky Variability Analysis (FAVA; Abdollahi et al. 
2017) and the 3FGL (Acero et al. 2015) monthly aper- 
ture photometry light curves! allow users to quickly ex- 
amine a source for relative flux increases, but do not 
provide flux calibrated characterization of a source. Do- 
ing so requires a full likelihood analysis of the region 
that takes into account the flux variations of all nearby 
sources, and generating a high cadence (e.g., daily) light 
curve using a full likelihood treatment over the entire 
duration of the mission can be very time consuming. 

To address this need, we developed the Fermi-LAT 
Light Curve Repository (LCR; Kocevski et al. 2021)?, 
consisting of a public database of light curves for vari- 
able Fermi-LAT sources on a variety of timescales. The 
repository provides publication-quality light curves on 
timescales of 3 days, 7 days (weekly), and 30 days 
(monthly) for 1525 sources deemed variable in the 
4FGL-DR2 catalog (Ballet et al. 2020). The repository 
consists of light curves generated through full likelihood 
analyses of the sources and surrounding regions, provid- 
ing flux and spectral index measurements for each time 
bin. Hosted at NASA's Fermi Science Support Center 
(FSSC), the LCR provides users with on-demand access 
to this light curve data, which is continually updated as 
new data becomes available. The repository is a new 
resource to the time-domain and multi-messenger com- 
munities for associating and monitoring LAT sources, 
in particular facilitating identification of time intervals 
with high fluence or flux. 


1 https: //fermi.gsfc.nasa.gov/ssc/data/access/lat /4yr_catalog/ 
ap_lcs.php 

2 https://fermi.gsfc.nasa.gov/ssc/data/access/lat/ 
LightCurveRepository / 
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Figure 1. Top: Histogram of the source types included in 
the LCR population, defined as having a variability index 
above 21.67, corresponds to «196 chance of steady emission. 
Center and bottom: The population of sources included in 
the LCR (empty markers). In grey we report the 4FGL-DR2 
population not included in the LCR as reference. The energy 
range of integration is 0.1—100 GeV. 


This paper is organized as follows. In Sec. 2 we 
describe the gamma-ray sources included in the LCR. 
Sec. 3 is devoted to the description and discussion of 
the automated data analysis employed by the reposi- 
tory. We discuss analysis caveats and data usage best 
practices in Sec. 4. We briefly summarize and conclude 
in Sec. 5. We provide a Quick Guide in Appendix A 
with the aim to help familiarize users with the LCR 
data portal. 


2. THE SOURCE SAMPLE 


Motivated by the science described above, the LCR 
focuses on the sources in the 4FGL-DR2 catalog that 
have variability indices greater than 21.67, where the 
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variability index can be thought of as a proxy for the av- 
erage fractional variability 6F/F, with ôF measured on 
timescales of 1 year. As defined in the 4FGL-DR2 cat- 
alog, which is based on 10 years of survey data, sources 
with such a variability index over 10 years are estimated 
to have a less than 1% chance of being steady. The re- 
sulting sample consists of 1525 sources, or roughly 26% 
of the 4FGL-DR2 catalog. The vast majority of these 
sources is blazars, further classified as flat spectrum ra- 
dio quasars (FSRQ), BL Lacs (BLL), and blazar candi- 
dates of unknown type (BCU), making up roughly 38%, 
31%, and 24% of the repository sample, respectively, or 
77%, 36% and 26% of their respective 4FGL-DR2 class. 
This is consistent with the fact that the LAT class with 
the largest number of variable sources are FSRQs. We 
also count 9 pulsars (PSR) and several sources of other 
kinds, as illustrated in the top panel of Fig. 1. In the 
middle and bottom panels of the same figure, we show 
plots of the spectral shape parameters (power-law index 
and total energy flux) and variability parameters (vari- 
ability index and fractional variability from DR2) for the 
sources included in the LCR (colored points) overplotted 
on the full 4FGL population (grey points). The former 
shows that most of the sources not included are fainter 
with harder spectra, i.e., mainly BL Lac objects. This is 
because, as can be seen in the bottom panel of the same 
figure, BL Lac objects are intrinsically less variable than 
FSRQs in the Fermi-LAT energy range. 


3. AUTOMATED DATA ANALYSIS 


Generating 3-day, 7-day, and 30-day light curves for 
each of these sources for over 13 years of data requires 
the analysis of over 3.7 million individual time bins. 
'The LCR analysis pipeline runs on a cluster hosted at 
SLAC National Accelerator Laboratory. Generating the 
light curves over the entire mission to date for the in- 
cluded sources requires approximately 3 months (7400 
core years). 

In the following subsections we present the details of 
the analysis and computational strategies used to gen- 
erate the LCR data products. 


3.1. Analysis Technique & Tools 


The characterization of LAT sources is typically per- 
formed using a maximum likelihood analysis (Abdo 
et al. 2009), in which the parameters of a model de- 
scribing the point sources and diffuse isotropic gamma- 
ray emission in a given region of the sky are jointly opti- 
mized to best describe the observed photon distribution. 
'The light curves of the LCR are obtained by performing 
an unbinned likelihood analysis, in which the full spatial 
and spectral information of each photon is used in the 
maximum likelihood optimization. 
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Figure 2. A model map, with 0.1 degree resolution, for 
a single weekly time bin of the region surrounding FSRQ 
4C +28.07. This source contains 15 other variable sources 
within a 12° radius, highlighting the need to model all vari- 
able sources within the ROI. 


The LCR analysis is performed with the standard LAT 
Fermitools? (version 1.0.5) using the P8R3_SOURCE_V2 
instrument response functions on P8R3_SOURCE class 
(Atwood et al. 2013; Bruel et al. 2018) photons selected 
over the energy range covering 100 MeV-100 GeV. Note 
that the energy dispersion is neglected in the unbinned 
analysis mode. This is not expected to impact the qual- 
ity of the analysis. However, repeating the analysis with 
a binned likelihood approach would return somewhat 
(mostly non significant) results. For each source and 
time bin, photons are selected from a circular region of 
interest (ROI) of radius 12? centered on the location of 
the target source. The ROIs are analyzed separately. By 
contrast, for the standard FGL catalogs, fluxes and spec- 
tral parameters of multiple sources are extracted from 
the same ROI. The size of the ROI is conservatively cho- 
sen to be as large as the 9596 containment radius of the 
LAT energy-dependent point-spread function (PSF) at 
100 MeV. Additional data selection cuts are imposed to 
exclude photons associated with regions and periods of 
known solar flares and gamma-ray bursts. A zenith an- 
gle limit of 90? to strongly reduce contamination from 
gamma rays produced through interactions of cosmic 
rays with Earth's atmosphere (Earth limb). 

Following the 4FGL-DR2 catalog, the LCR sources 
can have one of three different spectral types: 


3 https://fermi.gsfc.nasa.gov/ssc/data/analysis/software/ 


e Power-law (PL): 


dN/dE = No(E/Eo) * (1) 


e log-parabola (LP): 


dN/dE = Ny(E/E,)-(**Plos(5/E) (2) 


e and subexponentially cutoff power-law (PLEC): 


dN/dE = No(E/ Eg) ^ e "F^ (3) 


A majority of LCR sources are best described by a PL 
or LP spectral shape. The PLEC spectral shape best 
represents 11 additional sources, comprising the 9 pul- 
sars in the LCR sample and the two brightest FSRQs 
(CTA 102 and 3C 454.3). 

The normalization of the source spectrum in the 
model is left free to vary, while the spectral shape pa- 
rameters are initially fixed to their 4FGL-DR2 cata- 
log values. The model also includes all gamma-ray 
sources in the 4FGL-DR2 catalog within a radius of 
30? from the ROI center. The normalization of each 
variable source in the ROI is also left free to vary 
in the model, with spectral shape parameters fixed 
to their catalog values. In addition to the point 
sources, Galactic and isotropic background components 
are included in the model. The Galactic component, 
gll iem v07.fits!, is a spatial and spectral template 
that accounts for interstellar diffuse gamma-ray emis- 
sion from the Milky Way. The isotropic component, 
iso.P8R3.SOURCE. V3. v1, provides a spectral template to 
model the remaining isotropic events, including contri- 
butions from the residual charged-particle background 
and the isotropic celestial gamma-ray emission. The 
normalizations of both the Galactic and isotropic com- 
ponents are allowed to vary during the fit. The free pa- 
rameters of the model are varied to maximize the likeli- 
hood of observing the data given the model. An iterative 
fitting strategy, which varies the required fit tolerance? 
over three steps (1, 1074 and 1078), is employed to min- 
imize the number of time bins in which the likelihood fit 
does not successfully converge. Once fit convergence is 
achieved with the tightest tolerance, a second round of 


4 https: //fermi.gsfc.nasa.gov/ssc/data/access /lat / 
Background Models. html 


5 The fit tolerance is the relative convergence tolerance that is 
specified in the gtlike routine and passed to the optimization 
algorithm used to maximize the log likelihood function. The fit 
tolerance can be loosened to help achieve fit convergence. On the 
other hand, fits with tighter fit tolerances exhibit overall lower 
fractional errors on the resulting flux estimations. Therefore an 
iterative approach has been adopted for the LCR. 
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Figure 3. An example 3-day light curve as can be found in the LCR. It spans more than 11 years of LAT data and refers to 
the source at the center of Figure 2, FSRQ 4C +28.07. This light curve is obtained by specifying the analysis option that the 
spectral index of the source is free to vary. The data gap in 2018 is due to the temporary shut down of the instrument because 


of a solar panel anomaly. 


fitting is performed in which a spectral shape parameter 
of the target source is allowed to vary, namely, photon 
index (I) for PL, a for LP, and ^4 for the PLEC model. 
All other parameters remain fixed, i.e., 0 for LP, and a 
and y2 for the PLEC model. 

Figure 2 shows a model map for a single time bin 
of the region surrounding FSRQ 4C +28.07 (4FGL 
J0237.8--2848), which contains 15 other variable sources 
within a 12? radius. Figure 3 shows the resulting 3-day 
light curve spanning more than 11 years of LAT data 
for AFGL J0237.8+2848. 

A likelihood ratio test (Neyman & Pearson 1928) is 
used to quantify the significance of the target source 
above the background; specifically the test statistic 
(Mattox et al. 1996): 


TS = —2 log(Lo/ L). (4) 


The TS compares the maximum value of the likelihood 
function Lo evaluated for the parameter values that 
maximize the likelihood under a background-only null 
hypothesis (i.e., a model that does not include a tar- 
get source), with L, the likelihood function evaluated at 
the best-fit model parameters when including the target 
source. In the null hypothesis, TS is distributed ap- 
proximately as x? (Wilks 1938), and the analysis rejects 
the null hypothesis when the test statistic is greater than 
TS > 4, which is roughly equivalent to a 2 o rejection cri- 
terion for a single degree of freedom. Using this thresh- 
old value of TS as the detection criterion, the LCR em- 


6 The null hypothesis is tested against the presence of one source 


at a known position (at the center of the ROI) 


ploys a Bayesian profile likelihood method" to calculate 
the 9596 confidence level upper limits for any interval 
that yields a TS < 4, and also extracting a flux estima- 
tion of the target source for any interval with a TS > 1. 
The resulting best-fit values (or upper limits) for photon 
flux (cm? s^ 1), energy flux? (GeV cm"? s^ 1), and asso- 
ciated spectral shape are saved to the LCR database, for 
both the fixed and free photon index analyses. This pro- 
cedure ultimately allows for a user-selectable detection 
threshold and spectral fit method, as the flux estimates 
and upper limits between 1 < TS < 4 are both recorded 
in the LCR database, as well as the results from both 
the fixed and free spectral analyses. 

In Fig. 4 we compare, for the case of the bright quasar 
'Ton 599, the best-fit flux values found using a fixed spec- 
tral index vs. leaving the spectral index free to vary. The 
colors in the plot refer to the three cadences analysed. 
The results are generally stable, with greater uncertain- 
ties for lower fluxes for higher cadences due to lower 
statistics in the ROI in each time bin. 


3.2. Computational Strategy 


For analyses of a relatively small number (a few thou- 
sands) of photons, an unbinned likelihood analysis can 
be performed rapidly (a few minutes), but as the num- 
ber of events increases, the time to perform the analysis 
can become prohibitive. This limitation becomes in- 
creasingly burdensome when the need arises to perform 


T A description of the Bayesian profile likelihood method employed 
by the LCR can be found at https://fermi.gsfc.nasa.gov/ssc/ 
data/p7rep/analysis/scitools/python_tutorial.html. 


8 The energy flux is obtained from the photon differential flux dN 
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Figure 4. Comparison of the fluxes [ph/cm? /s] in time bins 
obtained with the fitting pipeline keeping the spectral index 
fixed versus letting the spectral index be free. Empty circles 
indicate results from analyses that did not converge, i.e., 
MINUIT Return Code q # 0. MINUIT (James 1994) is the 
optimizer adopted in the likelihood analysis. The plot shows 
the fluxes measured for the quasar Ton 599. In this case, the 
fit for the highest-flux point (with TS~ 4) converged when 
the spectral index was set free, but did not converge when 
the spectral index was fixed. 


a source characterization over a large number of time 
bins. A binned likelihood analysis could alleviate this 
issue. However, information is lost when binning data. 
'The LCR tackles the computational overhead by paral- 
lelizing the process of performing a full unbinned likeli- 
hood analysis. In order to produce a high cadence light 
curve over the entire lifetime of the mission in a reason- 
able amount of time, the LCR distributes the analyses of 
the light curve bins to separate nodes in a computer clus- 
ter hosted at the SLAC National Accelerator Laboratory 
and utilizing the IBM Spectrum LSF workload manage- 
ment platform. The parallelization allows for thousands 
of time bins to be analyzed simultaneously, with the net 
effect of drastically reducing the time to generate light 
curves over the entire duration of the mission. 


3.3. Exposure analysis 


The Fermi spacecraft has been executing a sky- 
scanning strategy for the LAT for the great majority of 
the mission, generally reaching an almost uniform full- 
sky coverage daily, in fact, every two orbits (ie., ap- 
proximately every three hours). During the mission the 
scanning strategy has been interrupted occasionally for 
targeted observations. The longest-duration such obser- 
vation was a modified observing strategy executed from 
2013 December to 2015 July that favored the Galactic 
center region. In March 2018, the seizing of the drive 
motor for one of the solar panels forced the temporary 


shut down of the LAT and a redefinition of the survey 
mode?. This caused the exposure to be limited to a part 
of the sky for some period of time. The effects of these 
periods of uneven exposure can be seen in the 3-day light 
curves. Therefore we quantify and discuss these effects 
in this section. 

The number of photons in a given ROI is modulated 
by the exposure of the observation, which is given by the 
product of the LAT livetime (during the brief readout 
time of a photon or cosmic-ray interaction the LAT, the 
instrument is ‘dead’ to triggering on other interactions) 
and the energy and angle-dependent effective area of 
the LAT. In general, analyses integrating over long ob- 
servation times can have small fractional exposure vari- 
ations across an ROI, but shorter-timescale analysis can 
be significantly impacted by low exposure for a particu- 
lar region of the sky. In fact, since the beginning of the 
post-anomaly modification of the sky-scanning strategy 
in February 2019, a number of the 3-day and 7-day ca- 
dence time bins have no exposure in particular portions 
of the sky. This is due to the constraints on the direc- 
tion of the zenith of spacecraft with respect to the Sun!9. 
Fig. 5 shows the range of values of the LAT exposure in 
the sky for each time bin of 3-day cadence. The maps 
in Fig. 6 illustrate the positions of all the sources ana- 
lyzed in the LCR, overlaid on the counts map (top) and 
an exposure map (bottom) for an typical 3-day cadence 
time bin with a zero minimum exposure: some sources 
fall in the region of the sky with zero exposure, which 
naturally translates into zero events to analyze in the 
ROI. 

In these cases, or in those with minimum exposure 
orders of magnitude lower than the maximum (shaded 
pink regions in the maps in Fig. 6), the pipeline typically 
returns an upper limit on the source flux. However, in 
some cases it could still find a solution that maximizes 
the likelihood, but warns of a poor fit quality or an un- 
reasonably low error value. For this reason, we recom- 
mend LCR users use caution when using the 3-day and 
7-day cadence light curve for these low-exposure inter- 
vals, or exclude the affected time bins entirely from their 
analyses. 

To illustrate the effect of low exposure on the photon 
counts, we generated binned all-sky counts and related 
exposure maps in the same time bins used for the LCR 
data analysis. We consider events in the 0.1-100 GeV 


9 https://fermi.gsfc.nasa.gov/ssc/observations/types/ 
post. anomaly/ 


10 Note that the ‘gaps’ in the exposure for different time bins are 


not always in the same portion of the sky 
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Figure 5. The gray band marks the range of LAT on-sky exposures for each time bin of 3-day cadence. The pink bands, in 
order, mark time ranges of the original rocking angle mission profile, the Galactic center monitoring, and the time gap between 
the solar panel anomaly and the beginning of the new survey mode. Note that some time bins have a minimum exposure value 
of zero: this means that some part of the sky was not observed during that time interval. The exposure maps are computed at 


the central energy in the 0.1-100 GeV band. 


(D Details on the Observatory sky-survey profiles can be found at https://fermi.gsfc.nasa.gov/ssc/observations/types/allsky/. 


energy range, and the maps are produced in HEALPix!! 
format with the tool gtexpcube2!?, using the same anal- 
ysis setup and IRFs as for the automated likelihood anal- 
ysis. The HEALPix format allows us to easily extract 
the average exposure in every ROI, and the pixel res- 
olution (HEALPix order 6) matches the resolution of 
the pre-computed livetime cubes provided for LAT data 
analyses. 

For each source considered in the LCR and each time 
bin of 3-day and 7-day cadences, we extract the average 
exposure and the total photon counts in the ROI cen- 
tered on the source between 0.1 and 100 GeV. In Fig. 7 
we illustrate the statistics of sources with a given num- 
ber of time bins in the light curve that had no photons 
in the ROI or had fewer than 20 photons (considered as 
an arbitrary threshold for low-statistics analysis). Hun- 
dreds of sources have up to 30 time bins with zero events 
in the light curve and/or less than 20 events, while tens 
of sources have fewer than 20 photons in more than 50 
time bins (note that the total number of time bins in 


11 http://healpix.jpl.nasa.gov/ (Gorski et al. 2005) 


12 The Fermitools, as part of the maximum likelihood calculation, 
automatically account for the exposure, as describet on the of- 
ficial mission web page at https:/ /fermi.gsfc.nasa.gov /ssc/data/ 
analysis/documentation/Cicerone/Cicerone_Data_Exploration/ 
livetime_and_exposure.html. The module gtexpmap is used to 
compute the exposure when performing unbinned analyses, as 
for the in the LCR analysis pipeline. In this section, we compute 
the exposure through the gtexpcube2 module, which is used 
when performing binned likelihood analyses. We stress that this 
is for purely illustrative purposes; the exposure maps computed 
in this section were not used in the LCR likelihood analyses. 


a 3-day cadence light curve is more than 1680 as of 
now). The monthly cadence, due to the longer inte- 
gration time, should not be affected by this issue; how- 
ever, we still recommend that users check the exposure 
of each bin used in their analyses. The exposures within 
the 12° radius ROIs for each source for the 3-day and 
7-day cadences are provided in the LCR downloadable 
data. 


4. USAGE CAVEATS 


In this section we provide a number of caveats to be 
mindful of when using the LCR data for scientific re- 
search. This list is also available and will be updated 
periodically on the LCR data portal. 


1. The LCR provides fit results from likelihood anal- 
yses that both did and did not converge. How- 
ever, it is important that the end user is aware 
that results from analyses that did not converge 
should be considered suspect and not be used for 
higher-level analyses (e.g., multi-frequency cross- 
correlation, power spectral density, or studies of 
variability). The convergence status for a partic- 
ular time bin is recorded in the MINUIT Return 
Code parameter and non-convergent analyses are 
hidden by default, but are optionally accessible to 
the user. A review of the current version of the 
LCR fit results for the first fourteen years of mis- 
sion data shows that the analyses for 0.796 of time 
bins, for all sources, did not converge when the 
spectral indices were held fixed, and that ~35% 
did not converge when the spectral indices were 
free to vary. 
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Figure 6. Top: Counts map for a 3-day time bin with 0 min- 
imum exposure; Bottom: Exposure map for the same time 
bin. The white pixels mark the positions of LCR sources. 
The light blue and magenta circles in the maps highlight the 
cases of two ROIs lying in a region of the sky with zero and 
non-flat exposures, respectively. 


2. From its definition, equation 4, the TS is a man- 
ifestly positive quantity. However, negative T'S 
values can sometimes be obtained when the pa- 
rameters reach the limits of their allowed intervals 
without having maximized the likelihood profile. 
Fit results obtained from intervals that resulted in 
negative TS values should be considered suspect 
and not used in higher-level analyses. The cur- 
rent version of the LCR data products for the first 
fourteen years of the mission for all LCR sources 
have only a few bins with negative TS results, per 
cadence, for both the fixed and free fits. 


3. While time intervals containing gamma-ray bursts 
(GRBs) and solar flares have been removed from 
the LAT data prior to the likelihood analyses, pos- 
sible contamination by the proximity of the qui- 
escent Sun has not been accounted for, nor have 
those time ranges been excluded. The angular sep- 
aration of the Sun from the target source is pro- 
vided for each time bin. A total of 175 GRBs and 
266 solar flare time intervals were excluded from 
the data prior to performing the maximum like- 
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|—— Weekly Counts < 20 
|— Weekly Counts = 0 
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Number of bins 
Figure 7. In light blue, we show the distribution of LCR 
sources by number of time bins with zero photons within 
the ROI used for their analysis. In magenta, we show the 
distribution of LCR sources according the number of time 
bins with fewer than 20 photons within the ROI used for 
their analysis. Darker shades of the lines refer to the 3-day 
cadence, while lighter shades refer to the weekly cadence. 
The low-exposure time bins represent the 0.196 of the total 
time bins. 


lihood analysis for the LCR first fourteen years 
of light curve data. This affects less than about 
0.0196 of the time bins for all the LCR sources per 
cadence. 


4. Because the LCR analyses are made available in 
real-time (new analysis results are generally made 
available within 24 hours of being processed), the 
results are not validated by the LAT Collaboration 
prior to release. Users are encouraged to perform 
sanity checks by examining the ratio of flux to flux 
uncertainty vs. the square root of the TS, and 
the distributions of fit results, e.g., flux, flux un- 
certainties, spectral indices mean values (photon 
index [ for PL, a for LP, and yı for the PLEC 
model) and their uncertainties. Some examples of 
these are shown in Fig. 8 for the specific case of 
the FSRQ 4C +28.07 (4FGL J0237.84-2848). The 
flux over the flux uncertainty ratio is expected to 
be approximately proportional to the square root 
of the TS. Any outliers should be either further 
investigated or removed before using the data for 
higher-level analyses, as should any extreme out- 
liers from the data distributions. 


5. The free-spectral-index light curves provided by 
the LCR were produced using a model of the per- 
tinent region of the sky for which only the spectral 
index of the target source is set free, and those of 
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Figure 8. Example of validation plots for the FSRQ 4C +28.07. The top panels show the case for which the spectral index 
(a for this source) is fixed. In all the plots we are considering the 0.1-100 GeV energy range. The bottom panels show the 
case for which the spectral index is free to vary. The middle panels show the distributions of flux statistical uncertainties. Flux 
distributions are shown in the right panels. In this example, distributions are good except for the outliers highlighted within 
green squares, which should be either further investigated or removed. In the histograms, dashed lines represent results from 
all time bins. Solid lines represent bins that did not converge. Notice that outliers are not necessarily the results from analyses 
that did not converge. 


all the other sources were fixed to the 4FGL-DR2 
catalog values. Therefore, contamination induced 
by possible changes in the spectral indices of the 
sources surrounding the target are not taken into 
account. For instance, sources undergoing bright 
flares have been seen to also experience dramatic 
changes in their spectral indices, including changes 
in the curvature, at the same time (e.g., harder- 
when-brighter behavior). Therefore, bright, vari- 
able sources in the ROI can induce this type of con- 
tamination and must be considered when a target 
source is in close proximity to any bright variable 
sources. 


6. Some erroneously high flux values during periods 


of zero or low exposure, often associated to a very 
small (or zero) error, have been found for several 
sources (see Section 3.3). We recommend check- 
ing the exposure value indicated for each time bin 
and source ROI in the provided tables. If the light 
curve has fluxes with very small fractional uncer- 
tainties, rather than flux upper limit, in the time 
bins with zero or low exposure, the reported fluxes 
and uncertainties should be considered unreliable 
and those time bins should be excluded from con- 
sideration. Time bins with zero error on the flux 
estimations are automatically hidden but are op- 
tionally available to the user. 


7. A bug in the make4FGLxml.py tool that the 


LCR uses to generate the ROI models was re- 
cently identified, in which the reported exten- 
sion of the extended sources with RadialGauss 
spatial profile is the 68% C.L. value instead of 
the sigma value used in the XML model. Thir- 
teen 4FGL-DR2 sources have this spatial profile, 
namely, Crab IC, IC 443, Monoceros, HESS J1303- 
631, FHES J1501.0—6310, FHES J1626.9—2431, 
FHES J1723.5—0501, HESS J1825—137, W A1, 
Cygnus Cocoon, FHES J2129.9+5833, FHES 
J2208.4+6443 and FHES J2304.0+5406. None of 
these sources is deemed variable in the 4FGL-DR2 
catalog. However, they are present in the ROIs of 
101 LCR sources, which might result in a system- 
atic bias in the target flux values across the whole 
light curve. Correcting this issue will require re- 
processing the data for these sources, which is al- 
ready underway. Meanwhile users should be mind- 
ful of this issue when considering sources in the 
vicinities of these extended sources. 


. Particular care should be taken when using the 


light curve for the synchrotron component of the 
Crab. The Crab has different components treated 
as separate entries in the 4FGL-DR2 catalog: the 
extended emission, the inverse Compton emission 
from the PWN and the synchrotron emission from 
the pulsar. The LCR analysis generates only the 
light curve for the synchrotron (variable) emission, 
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while keeping the parameters of the other two com- 
ponents frozen in the fit. This could result in some 
contamination of the synchrotron component light 
curve deriving from any unaccounted-for variabil- 
ity of the other Crab components. On the other 
hand, the LCR approach overrides the apparent 
(not real) variability reported for the pulsar in 
4FGL-DR2 (and all other versions of the 4FGL 
catalog). 


5. PROSPECTS AND CONCLUSIONS 


The development of the Fermi-LAT LCR was moti- 
vated by the need for a coherent collection of light curves 
of variable gamma-ray sources observed by the Fermi- 
LAT in support of the time-domain and multi-messenger 
communities. By continuously reporting the flux evolu- 
tion and transition to high-flux states for many variable 
sources, the LCR is a valuable resource for triggering ob- 
servations of other observatories. Furthermore, the LCR 
can be used to validate the study of variable activity in 
neighboring faint sources, helping to identify potential 
contamination from flaring activity of a bright source. 
In this manuscript we described the automated analysis 
pipeline which will continuously update the repository 
with new data as soon as new observations by the LAT 
become available. We invite the community to use the 
LCR data products, and report any issues or suggestions 
to the LCR contacts at the Fermi Science Support Cen- 
ter. A brief guide for navigating the LCR to navigate 
the web site is provided in Appendix. 

In the future, the LCR source list and the resulting 
light curve data will be updated with every new Fermi- 
LAT source catalog, with the next LCR version planned 
to be released alongside the 5FGL catalog. Due to the 
computational expense of re-analyzing the full-mission 
light curves for the entire LCR sample, the LCR sample 
will not be updated for each catalog sub-release (e.g., the 
newly available 4FGL-DR3 (Abdollahi et al. 2022b) and 
the upcoming 4FGL-DR4). Currently, the LCR does not 
provide permanent identifiers that allow distinguishing 
between from different versions due to data reprocess- 
ing. This will be implemented in a future version of the 
repository. 


In conclusion, in this era of large surveys, the Fermi- 
LAT is the only high-energy gamma-ray observatory 
to continuously monitor variable sources, providing the 
all-sky coverage to identify gamma-ray counterparts to 
transient events at other wavelengths. We expect that 
the LCR will greatly enhance the usefulness of LAT 
data to the time-domain, multi-messenger and multi- 
wavelength communities. 
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APPENDIX 


A. QUICK USER GUIDE 
In this appendix we provide a comprehensive list of the main features of the LCR website!3. 


e The main page of the website features an interactive Catalog Map plotting the positions of all 4AFGL-DR2 
sources. Optionally, additional data may be overlaid on the map , e.g., real time Sun or Moon position, positions 
of IceCube neutrino alerts, and GRB error circles as reported in the Second LAT GRB catalog (Ajello et al. 
2019). The 1525 LCR sources are highlighted in dark gray, while by default the non-variable 4FGL-DR2 sources 
are marked in light gray. Hovering over any source displays a tooltip box showing its name and key characteristics 
as well as linking to its 4FGL light curve and spectrum, related FAVA entry, and LCR light curve if applicable!*. 
Below the map, a table is shown listing the 4FGL sources and important parameters. Clicking the name of a 
source included in the LCR opens a separate page dedicated to that source. 


e A Map Options menu provides numerous options related to the display of the Catalog Map. Options are provided 
to change the coordinate system and celestial projection. The marker label and color can be changed, as can the 
meaning of its size to indicate the variability index, average significance, or time-resolved significance in 3-day 
bins. 


e A Catalog Search toolbox allows the user to search for a specific source by name or Right Ascension and 
Declination. The results of the search are highlighted in the Catalog Map. Clicking on the linked name of 
the target source opens a dedicated page for the source in a new tab. 


e A Data Overlays toolbox allows the user to visualize a number of additional catalog overlays in the sky map. 
These catalogs include: the Fermi-LAT Gamma-ray Burst Catalog (2FLGC; Ajello et al. 2019), the IceCube 
Neutrino Alerts!?, and the FAVA Flare Catalog (2FAV; Abdollahi et al. 2017). Activating any of the catalogs 
will also add the related table under the map. 


e The dedicated page for each source shows an interactive light curve displaying the detections and upper limits for 
that source. Options are provided to show the 3 day, 7 day (1 week), or 30 day (1 month) cadence light curve. By 
default only the significant flux points (with TS>4) are shown, while upper limits are shown for less-significant 
time bins. However the user can choose to change the minimum detection significance through the drop-down 
menu, by selecting among the available options (TS,,;,—4,3,2,1). The Spectral Fitting option allows the user to 
choose to visualize either the best-fit values obtained with the fixed spectral index fit or the ones resulting from 
the fit with spectral index free to vary. A table listing the main characteristics of the selected source is provided. 
Additional information about the fit convergence, fit tolerance and detection ratio are reported for diagnostics 
purposes below the light curve plot. 


e A Data Download toolbox on each dedicated source page provides all data for that source for download. The 
data are provided in CSV and JSON formats!?. All data points are provided, potentially including unconstrained 
and possibly TS« 0 data points, or from analyses that did not converge. Guidelines for cleaning the data before 
use in analysis are given in Sec. 4. 


e Finally, the LCR contains a Usage Notes page which reviews the data analysis and modeling details, fitting 
strategy, and caveats for usage (see Sec. 4). 


13 https: //fermi.gsfc.nasa.gov/ssc/data/access/lat / 
LightCurveRepository/ 

14 This option can be disabled by unchecking the Source info option 
from the Map Options menu. 

15 https: //gcn.gsfc.nasa.gov/amon.html 

16 A description of the file formats can be found at https:// 
fermi.gsfc.nasa.gov/ssc/data/access/lat /LightCurveRepository / 
table. description.html 


